Stripe-hexagon competition in forced pattern forming systems with broken up-down 

symmetry 
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We investigate the response of two-dimensional pattern forming systems with a broken up-down 
symmetry, such as chemical reactions, to spatially resonant forcing and propose related experiments. 
The nonlinear behavior immediately above threshold is analyzed in terms of amplitude equations 
suggested for a 1 : 2 and 1 : 1 ratio between the wavelength of the spatial periodic forcing and the 
wavelength of the pattern of the respective system. Both sets of coupled amplitude equations are 
derived by a perturbative method from the Lengyel-Epstein model describing a chemical reaction 
showing Turing patterns, which gives us the opportunity to relate the generic response scenarios to 
a specific pattern forming system. The nonlinear competition between stripe patterns and distorted 
hexagons is explored and their range of existence, stability and coexistence is determined. Whereas 
without modulations hexagonal patterns are always preferred near onset of pattern formation, single 
mode solutions (stripes) are favored close to threshold for modulation amplitudes beyond some 
critical value. Hence distorted hexagons only occur in a finite range of the control parameter 
and their interval of existence shrinks to zero with increasing values of the modulation amplitude. 
Furthermore depending on the modulation amplitude the transition between stripes and distorted 
hexagons is either sub- or supercritical. 

PACS numbers: 82.40.Ck, 47.20.Ky, 47.54.+r 



I. INTRODUCTION 

External periodic forcing provides a powerful tool to 
analyze the response behavior of nonlinear pattern form- 
ing systems, allowing for instance the study of their inher- 
ently nonlinear mechanism of self organization. The early 
investigations of effects dealing with spatially periodic 
forcing were devoted to pattern formation in hydrody- 
namic systems followed by chemical syste ms [II [4 bl 01 La 
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The situation is comparable for temporal forcing of pat- 
terns whereas re- 
cent investigations on spatiotemporal forcing have mainly 
been motivated by chemical reactions [33J, |3J, |2a, 36] . 

In thermal convection, when deviations from the 
Boussinesq approximation come into play |3?1 l3S| . as 
well as in several chemical reactions of the Turing type 
[3^. lioj , the up-down symmetry of the fields describing 
the patterns, u(r,t) — > — u(r, t), is broken. In such two- 
dimensional extended systems rotationally symmetric in 
a plane, hexagons are the generic pattern close to thresh- 
old. In systems with a strong axial anisotropy, such as in 
nematic liquid crystals, this rotational symmetry is bro- 
ken [4ll lia . liM and hexagonal convection patterns 
do not occur. Thus the interplay between different bro- 
ken symmetries, such as the broken up-down symmetry 
with a weak anisotropy, leads to an interesting competi- 
tion between stripe and hexagonal patterns as has been 
shown recently 0, EE ^3 . 

Which kind of effects may be expected if besides the 
up-down symmetry the translational symmetry is simul- 
taneously broken by periodic modulation in one spatial 
direction? As the direction of the modulation wave- 



number defines a preferred direction the system becomes 
anisotropic. In the following we present the symme- 
try adapted amplitude equations and analyze the related 
hexagon-stripe competition in terms of these equations. 
In thermal convection this symmetry breaking may be 
achieved by modulating the container boundaries or the 
applied temperature difference periodically in one direc- 
tion 0,1a, UM- In chemical reactions however, the forcing 
is introduced by a spatially modulated illumination of 
the system [2E E3- The specific manifestation of the 
modulation-induced symmetry breaking in generic am- 
plitude equations depends on the ratio between the mod- 
ulation wave number k m and the wave number k c of the 
respective pattern. We will study two cases of resonant 
modulation: k m = 2k c and k m — k c . 

Exemplarily we will focus in this work on the Lengyel- 
Epstein model for a chemical reaction. This model, in- 
troduced in Sec. [HJ displays in some parameter range a 
Turing instability of the homogeneous basic state as will 
be shown in Sec. lllll Immediately above the threshold of 
a supercritical bifurcation the amplitudes of the unsta- 
ble pattern forming modes are still small enough to apply 
the powerful perturbational technique of amplitude equa- 
tions. In Sec. II VI we derive the amplitude equations from 
the given Lengyel-Epstein model in the limit of small 
spatially periodic modulations of the illumination acting 
on the chemical reactions. We consider a resonant ra- 
tio of 1 : 1 as well as of 1 : 2 between the wavenumbcr 
of the Turing instability and the modulation wavenum- 
bcr. The amplitude equations obtained in each of these 
resonant cases are generic and may be found for other 
systems having the same symmetry properties as the con- 
sidered Lengyel-Epstein model, for instance modulated 
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convection systems. Without spatially periodic forcing 
it is known that the resulting amplitude equations favor 
hexagonal patterns close to threshold. In Sec. we in- 
vestigate how the spatial forcing modifies the transition 
scenario between stripe and hexagon-like patterns and 
how the initial hexagonal patterns get distorted. Analyt- 
ical calculations are confirmed by numerical simulations 
in Sec. I VII Concluding remarks and an outlook on trav- 
eling stripe forcing within the scope of the underlying 
chemical reaction are given in Sec. IV 111 



II. THE LENGYEL-EPSTEIN MODEL 

As a basic model of a chemical reaction we choose the 
Lengyel-Epstein model as described in Ref. |4^. If the 
effect of illumination on the chemical reaction is taken 
into account, as done in Ref. |5Cj . the model may be 
described by the following equations 



dtu = a — cu — 4 
dtV = a | < '/ 
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In these equations the dimensionless concentrations of 
two relevant chemical substances are labelled with u(r, t) 
and v(r, t), whereas a, c, a and d denote dimensionless pa- 
rameters of the chemical system. Finally the effect of ex- 
ternal illumination is described by the field <I>(r, i), which 
can be identified as the control parameter of the system 
and decomposed into a spatially homogeneous and a spa- 
tially periodic part HJ El H3 



$(r,<) =4> + M{r) 



(2) 



For the second term of Eq. J2J) we assume a periodic 
dependence in one spatial direction as described by 



M(r) = 2Gcos(fc m x) 



(3) 



For a spatially homogeneous illumination Eqs. have 
a stationary and homogeneous solution given by 
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III. THRESHOLD FOR PATTERN FORMATION 

In a certain parameter range the spatially homoge- 
neous state described by Eqs. |@J becomes unstable 
against infinitesimal perturbations ui(r, t) and vi(r,t). 
In order to determine this interval the basic state is sep- 
arated from the inhomogeneous contributions to both 
fields by the ansatz 



u(r,t) 
v(r, t) 



m + ui(r, t) , 
v + vi(r,t) . 



(5a) 
(5b) 



In a next step the linerization of the basic Eqs. with 
respect to those small contributions ui(r, t) and vi{v,t) 
yields linear equations with constant coefficients, which 
may be solved by the mode solutions 



Mi(r, t) = u\e 



AH 



vi{v,t) = vie 



Af+ik-r 



(6) 



and thus leads to a set of homogeneous linear equations 
for the amplitudes ui and vi. From these equations 
the dispersion relation X(k) may be specified. It turns 
out that a strong and homogeneous illumination of the 
chemical reaction suppresses inhomogeneous perturba- 
tions. This is illustrated by Fig.^whercin the dispersion 
relation A(fc) is plotted for three different cases of the 
illumination strength 4>. 
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FIG. 1: The dispersion relation X(k) for a supercritical value 
4> < <f> c (dotted line) of the illumination strength cf>, for a 
critical value <j> = cj> c (solid line) and for a subcritical value 
4> > 4> c (dashed line). Parameters are a — 16, c = 1.0, d — 
1.5, o = 301. 



From the neutral stability condition 
Re(A(fc)) =0, =S> 0o (fc) 



(7) 



the neutral curve (f>o{k) may be determined, which has 
to be done numerically in our case. Some examples of 
characteristic neutral curves are plotted in Fig. [21 for dif- 
ferent values of the diffusion coefficient d. For d = 1.2 
the homogeneous state is either stable for a small or for a 
strong illumination <fi and is thus unstable only within a 
small island enclosed by the corresponding neutral curve 
in the </>-fc-plane. However, for d > 1A the homogeneous 
state remains only stable for strong illuminations. In this 
case the neutral curve <po(k) takes its maximum at the 
wavenumber 



1 /4«b(ug- 1) 



14 



(1 



,2^2 



(8) 



and the value of the illumination at the maximum of the 
neutral curve defines the critical illumination strength 
4>c = 4>o{k c )- This critical value is plotted in Fig. [21 as a 
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FIG. 2: The neutral curve <po(k) is determined numerically via 
the neutral stability condition in Eq. for different values 
of the diffusion parameter d — 1.2, 1.3, 1.6 (from inside to 
outside). The remaining parameters are a — 16, c = 1.0, <r = 
301 in all cases. The homogeneous state is unstable inside of 
each curve. The dotted line shows k c as a function of d. 



IV. AMPLITUDE EQUATIONS AND LINEAR 
STABILITY OF THEIR SOLUTIONS 

The amplitudes u\ and v\ of the linear modes described 
by Eq. © as well as the amplitudes of any other pattern 
are restricted by the nonlinear terms proportional to u% 
and vi in Eqs. Q). In two-dimensional systems with a 
broken up-down symmetry the common patterns close to 
threshold are stripes, squares or hexagons [&2]- Immedi- 
ately beyond a supercritical bifurcation, where the am- 
plitude of an emerging pattern is still small, the concept 
of the so-called amplitude equations is a very successful 
one to characterize the nonlinear behavior of patterns, as 
it is exemplified for several physical, chemical and biolog- 
ical systems in Refs. |52l l53|. The amplitude equations 
are obtained by a perturbative method from the basic 
Eqs. 0J, whereby the coefficients occurring in these am- 
plitude equations reflect the dependence of the patterns 
on the parameters of the specific system. 

A. Scheme of the derivation of amplitude equations 



function of the diffusion coefficient d and for various val- 
ues of the parameter c. The occurring Turing instability 
vanishes, when the inner closed curve in Fig. [21 shrinks to 
zero, i.e. below d c = 1.179 for the parameter set chosen 
in this figure. Below this critical value d c the station- 
ary homogeneous state becomes stable with respect to 
an oscillatory instability. 
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The small parameter used in the perturbative deriva- 
tion of the amplitude equations is the relative distance 
to the threshold <b r 



(9) 



Since the up-down symmetry for inhomogeneous field 
contributions is broken, i.e. U\,Vi -/-* —U\,—Vi, we 
expect hexagonal patterns to occur in some parameter 
ranges. Close to threshold hexagonal patterns are de- 
scribed by a superposition of three waves (stripe pat- 
terns) with wave vectors enclosing an angle of about 120° 
with respect to each other. If the two fields u± and v\ are 
rewritten in terms of a vector Wi = (ui,vi) T , the gen- 
erating field of hexagonal patterns may be represented 
by 



Wi 



£ q ) (A ie lkir + A 2 e^ r + A 3 e l ^ r ) + c.c. , 



(10) 



where the wave vectors (i = 1, 2, 3) of the three un- 
derlying stripe patterns read 



'" ' J) and k2 ' 3 = y (±73 



(ii) 



FIG. 3: The critical illumination (j> c = (j>o{k c ) is plotted as a 
function of the diffusion coefficient d for different values of the 
parameter c, with c = 0.6, 0.7, 0.8, 1.4 from top to bottom 
and a — 16, a = 301. The dotted line corresponds to that in 
Fig. El 



and the second component Eq of the eigenvector is given 
by 



E = -(c + kl) 1 -±^-^ 1 -^4 
v cj 4w u Q 1 + ul 



(12) 



The coupled amplitude equations deduced for the three 
envelope functions Ai (i = 1, 2, 3) are identical apart from 
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the term evoked by the forcing which naturally differs 
whether we choose in Eq. © a ratio 1 : 2 = k c : k m 
between the modulation wave number k c and the critical 
wave number k m or a ratio 1 : 1 = k c : k m . In the 
first case the modulation amplitude G should vary as 
G = G\ oc e, in the latter one as G = G2 oc e 3 / 2 . 
Choosing the ansatz specified in Eq. 1|10J) and expanding 
Eq. Q with respect to powers of the small parameter e as 
well as the small amplitudes Ai one obtains the following 
set of nonlinear equations for the three amplitudes Ai 



T d t A 1 =sA 1 + 6 A* A* + GxA\ + G 2 

~( ri \A l f+p\A 2 \' i +p\A^)A l: (13a) 
T Q d t A 2 =eA 2 + 8AlA\ 

-( 1 \A 2 \ 2 + p\A 3 \ 2 + p\A 1 \ 2 )A 2 , (13b) 
T d t A 3 =eA z + 8A\Al 

-( 1 \A 3 \ 2 + p\A 1 \ 2 + p\A 2 \ 2 )A 3 , (13c) 



which are quite similar to those presented in Refs. 35, 45, 
l54) | . Depending on which of the resonant cases mentioned 
above is considered either G\ or G 2 is nonzero. For fur- 
ther details on the derivation scheme of the amplitude 
equations we refer to Ref. |55| . where likewise reaction 
diffusion equations were studied in detail. 
The analytical expressions of the coefficients To, 8, 7, p 
appearing in the amplitude equations (|13|l as functions 
of the parameters of the basic equations Q are rather 
lengthy. Instead of giving their analytical forms we there- 
fore plot them in Fig. 0] as functions of the diffusion coef- 
ficient d while assuming different values for c. Rescaling 



as well as the coefficients of Eqs. (|13fl as follows 
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(14) 



and yields a simpler form of Eqs. (|13|l 

d- t A x = + SA* 2 A* 3 + Gi A{ + G 2 

- (\A X \ 2 + p\A 2 \ 2 + p\A z \ 2 ) A x 



d- t A 2 = v A 2 + M* 3 Al 



(\A 2 \ 2 + p\A 3 \ 2 ■ 



(\A 3 \ 2 + p\A 1 \ 2 ^ 



■p\Ai\ 2 )A~2: 



P\A 2 \ 2 ) A 3 



(15a) 
(15b) 
(15c) 



with 8 = ±1. Choosing tq = 7 = 1 and 8 = ±1 in 
Eqs. I|13l) reduces them to the same form as Eqs. <|15[) 
which besides the modulation amplitude Gi (i = 1,2) 
and the control parameter rj only depend on the nonlinear 
coefficient p. In most cases this coefficient takes values 
in the interval p S [1,2] as shown in Fig. [SI where p 
is plotted as a function of the diffusion coefficient d for 
different values of c. 




FIG. 5: The rescaled coefficient p = p/7 as a function of d 
for various values of c, with c = 0.7,0.9, 1.1, 1.3 from left to 
right. 



FIG. 4: The coefficients of the amplitude equations in 
Eq. 1131 are shown close to the 1 : 2 resonance, i.e. Gi = 0, 
as functions of the diffusion coefficient d for different values 
of the parameter c, from left to right c = 0.7, 0.9, 1.1, 1.3. 

of the amplitudes Ai (i = 1, 2, 3) allows to express time 



B. Functional 

The coupled nonlinear equations l|13|) may be consid- 
ered as a potential dynamics 



t d t Ai 



8T 
"8A* 



(16) 
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with the functional 

T=j2{l\A^-e\A^) + p -j^\A^\A^ 

-^-(A 2 + A* 1 2 )-G 2 (Ai+Al) 

- S (A 1 A 2 A 3 + A* 1 A* 2 A* 3 ) . (17) 

This functional or its rescaled version yielding Eqs. 115(1 
provides a useful criterion to decide which kind of pat- 
tern is most likely to occur in a given parameter range. 
Generally the pattern leading to the lowest value of the 
functional is preferred. Once again either G\ or G 2 is 
nonzero, corresponding to one of the considered resonant 
cases. 



one considers the stability of stripes in a 2 : 1 resonant 
case, the amplitude B — and the marginal stability 
condition leads to a polynomial in rj 

(P-1)V 2 + [2pQi (p-1)-^ 

-Gi (5 2 -p 2 Gi) =0, (21) 

whose solutions are 

2pGi (1 - p) + S 2 
** = 2(l_-0» 

± * ^ v^+^a-p). (22) 

Consequently single mode solutions (stripes) are linear 
unstable within the region bounded by those values of r\ 
and linear stable in the complementary area. 



C. Linear stability of hexagons and rolls 

Since we only consider a spatially periodic modulation 
in x-direction, the translational invariance in y-direction 
is preserved. This holds for both resonant cases under 
investigation. Therefore two out of three amplitudes, 
namely A 2 and A3 are degenerated and the remaining 
two stationary amplitudes, called A and B, are deter- 
mined by two of the three Eqs. I|15fl . The linear stabil- 
ity of the stationary and nonlinear solutions is analyzed 
by superimposing small perturbations to the latter ones, 
meaning we use the following ansatz: 



A 1 =A + f x 1 (t), A 2 , 3 = B + p 2 , 3 (t) 



(18) 



Accordingly the amplitude equations (|15f) are linearized 
with respect to the small perturbations pi (i = 1,2,3) 
which leads to the following set of equations 

%n = (77 - 2\A\ 2 - 2p\B\ 2 ) m + (Gi - A 2 ) pi 

~pAB*(p 2 +p 3 )+ (SB* - pAB) (,4 + /J3) , 

(19a) 

d- t p 2 = ( v - (2 + p)\B\ 2 - p\A\ 2 ) p 2 - B 2 p* 2 

- pA*Bp x + (SB* - pAB) p{ 

- p\B\ 2 p 3 + (SA* - pB 2 ) pi , (19b) 
%i 3 = (ry - (2 + p)\B\ 2 - p\A\ 2 ) p 3 - B 2 p\ 

- pA*Bp l + (SB* - pAB) p{ 

- p\B\ 2 p 2 + (SA* - pB 2 ) f4 . (19c) 

These coupled equations may be transformed into a set 
of homogeneous linear equations by means of an ansatz 
for the perturbations 



at 1 * cr* t 

t* = v u e + v 2i e 



(» = 1,2,3) 



(20) 



which allows to determine the growth rate a(k, . . .) 
of those perturbations. Marginal stability, i.e. 
Re[er(fc, ...)]= 0, then yields the stability boundary as a 
function of the control parameter r\. As an example, if 



V. NONLINEAR SOLUTIONS OF THE 
MODULATED AMPLITUDE EQUATIONS 

In this section the bifurcation-scenarios from the ho- 
mogeneous basic state and thus the weakly nonlinear 
solutions of Eqs. 1 (13) 1 are discussed for three different 
cases: first of all we consider an unmodulated system 
(Gi = G 2 = 0), then we study the modulated cases al- 
ready taken into account in Sec. lIVI namely either G\ 7^ 
and G 2 = or Gi = and G 2 ^ 0. Without modula- 
tions one has either a single mode solution with modu- 
lus \A\\ — A and \A 2 \ — \A 3 \ = 0, or hexagons corre- 
sponding to a three mode solution with coinciding mod- 
uli I Ai| = IA2I = I ^3 1- In ca se of a finite modulation 
amplitude single mode solutions as well as three mode 
solutions with amplitudes satisfying |Ai| 7^ IA2I = IA3I 
may be encountered. By means of the determined nonlin- 
ear solutions the functional value of the latter ones may 
be calculated through Eq. I|17|) . which gives rise to the 
intervals wherein either the one mode or the three mode 
solution is preferred. 



A. Solutions of the unforced amplitude equations 



Without modulations Eqs. (|15fl have cither a single 
mode solution |A| = ^Jrj or a three mode solution with 
equal moduli, i.e. with amplitudes fulfilling |Ai| = 
IA2I = I A3 1 = A. Assuming a relative phase shift 9j for 
each linear mode contributing to the three mode solution 



Aj = Ae l9) 



(23) 



one obtains a nonlinear equation for the common modu- 
lus A 

= 7]A + 5e lB A 2 - (1 + 2p)A 3 , (24) 

with 8 = ±1 and the sum of the phase angles 6 = B\ + 
6> 2 + #3 • Eq. I|24l) has two real solutions 

1 



A+ = 



2(1 + 2p) 



S ±J S 2 +4 V (l + 2p) 



(25) 
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where A + corresponds to the larger one of the two am- 
plitudes with 8=1, while A_ is obtained for 8 = — 1. If 
8 = 1 the phase angle is 9 = and one expects regular 
hexagons whereas 5 — — 1 yields a phase angle 9 = 7r 
corresponding to inverse hexagons. Comparing the func- 
tional given in Eq. (|17|l of those solutions one may as- 
certain: In case 6=1 regular hexagons are preferred as 
Tjj < TT, and if 8 = — 1 inverse hexagons are most likely 
to occur because their functional value is the lowest one, 
i.e. < T~^. 

B. Solutions of the forced amplitude equations for 

k m = 2fc c and Gi / 0, G 2 = 0. 

In the case of a 2 : 1 resonance with a finite modula- 
tion amplitude G\ two of the three amplitudes coincide 
A 2j3 = Be 1023 , the third one being A 1 = Ae m , Thus 
the remaining two amplitudes A and B are determined 
by the following two nonlinear equations 

(rj - A 2 - 2pB 2 ) A + G\A + 8e iB B 2 = , (26a) 
(77 - (1 + p)B 2 - pA 2 ) B + 8e m AB = , (26b) 

with 9 = 81 +02 +03- The single mode solution generated 
by those coupled equations is given by 

A=y/r l + G 1 , B = (27) 

and bifurcates supercritcally at the threshold e = — G\ 
or rj = — G± in rescaled units. A linear stability analysis 
of the single mode solution according to the scheme de- 
scribed in Sec. II V 01 shows that the single mode is linear 
unstable in the region delimited by the two values of rj 

2pGi (1 - p) + S 2 

m ' 2 = 2(1-0' 

± * , 2 ^ 2 +4G 1 (l-p). (28) 

The three mode solution with 8 e l ® = 1 is generally 
preferred and the appendant amplitudes A and B may 
be determined by solving Eqs. (j2*FJ|> with a standard com- 
puter algebra program. Instead of presenting the rather 
complex analytical solutions we opt for the characteristic 
bifurcation diagrams based on four different parameter 
sets in Fig. EI which ought to cover the typical bifurcation 
scenarios encountered. In Fig. [5] the linear stable inter- 
vals of each branch are visualized by solid lines, whereas 
unstable ones are represented by dashed lines. 

Part (a) of Fig.E]shows the well known bifurcation sce- 
nario of the unmodulated case (G\ = 0) for supercritcally 
bifurcating stripes and subcritically bifurcating hexago- 
nal solutions [54J. With increasing values of the modu- 
lation amplitude G\ the degeneracy of the three mode 
solution vanishes and only the onset of the single mode 
solution starts at rj c = — G%. This removal of degeneracy 
is already best reflected in Fig.EJb), wherein G\ = 0.05: 




-0.11 -0.09 fj -0.30 0.00 0.30 



FIG. 6: The amplitude of supercritically bifurcating stripes 
and the subcritically bifurcating hexagons with three degen- 
erated amplitudes are shown in (a) as a function of n for the 
unmodulated case. In each case the solid lines mark the linear 
stable part of each branch and the dashed parts the unstable 
ones. For finite modulation amplitudes Gi the three mode so- 
lution with A always corresponding to the upper branch and 
B to the lower branch, as well as the amplitude of a single 
mode solution are plotted as a function of the control param- 
eter n for three different values of the modulation amplitude 
Gi = 0.05 (b), Gi = 0.11 (c) and Gi = 0.25 (d) whereas the 
nonlinear coefficient p is fixed at p = 1.8 and 6 = 1. Each 
vertical dotted line reflects the n- value where the values of the 
functional T of the one and the three mode solution coincide. 



the subcritical branches of the hexagonal solution, A cor- 
responding to the larger amplitude and B to the smaller 
one, do not overlap anymore. For finite values of G\ the 
upper hexagon-branch, i.e. the A-branch, always emerges 
from the single mode solution, while the B-branch which 
is the lower hexagon-branch bifurcates from zero with de- 
creasing values rj. If G\ is small enough the bifurcation 
of the three-mode solution is subcritical as depicted in 
Fig. [S] (b) and Fig. HJc), whereas for large modulation 
amplitudes this bifurcation is shifted to smaller ?7-values 
while becoming supercritical as sketched in Fig.El(d). 

It turns out that for intermediate values of G\, such 
as those used in Fig.[|Jb) and Fig.^c), the three mode 
solution bifurcates subcritically from the single mode so- 
lution not only for small values of the control parameter 
•q but also for large values of 77 as it may be concluded 
from Fig. For large values of G\, e.g. G\ = 0.25 in 
Fig. Eld) , the lower bifurcation point becomes supercrit- 
ical while the upper one nevertheless remains subcritical 
at first. By comparing the bifurcation diagrams shown 
in Fig. [7| and Fig. HJd) it can easily be seen that the 
existence range of the three mode solution shrinks with 
increasing values of the modulation amplitude. Simul- 
taneously the hysteresis becomes less pronounced with 
increasing values of G\ changing the subcritical bifurca- 
tion to a supercritical one at small 77-values as well as 
the upper bifurcation point from subcritical to supercrit- 
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FIG. 7: The amplitudes A(rj) and -B(??) are plotted for Gi — 
0.05 and Gi = 0.11 (already shown in Fig. |SJ over a broader 
interval of rj which reveals that the upper bifurcation point 
of the three mode solution moves for increasing modulation 
amplitudes Gi to decreasing values of rj. As in Fig. |S|p = 1.8 
and 5 = 1. 

ical for even larger values of G\. This latter point can 
be nicely seen in Fig. [S] which shows that beyond some 
critical value Gi c of the modulation amplitude the three 
mode solution ceases to exist. Since both bifurcations 
from the single to the three mode solution are already 
supercritical before their domain of existence shrinks to 
zero, the critical value G\ c may be determined using the 
condition that the two ri-values given by Eq. I|28|l coin- 
cide, i.e. Tji = rj2- Hence for modulation amplitudes G\ 
greater than the critical value 



single mode solutions are stable for all rj and three mode 
solutions do not exist any more. Since Fig.[S]evinces that 
in most cases p takes values larger than one, G\ c should 
always be a positive parameter. 

The domain of existence of three mode solutions, the 
stability range of single mode solutions as well as the 
points of coincidence of the corresponding functionals 
(vertical dotted lines) are represented in Fig. Etb)-(d) 
and Fig. for three different values of G\, These rel- 
evant properties are recapitulated as a function of the 
modulation amplitude G\ in Fig. |S| Within the region 
bounded by the dotted lines three mode solutions ex- 
ist. Along the dashed line the functional values of one 
and three mode solutions match and the area lying in 
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FIG. 8: Between the upper and lower dotted lines (the lower 
one overlapping partially with the solid line) lies the domain 
of existence of three mode solutions. The functionals of sin- 
gle and three mode solutions coincide along the dashed line 
(which overlaps once again partially on its lower branch with 
the solid one) and the region lying in between these lines is 
the one where three mode solutions are favored. Single mode 
solutions are linear unstable inside the area delimited by the 
solid line whereas linear stable outside and the dashed-dotted 
line corresponds to the threshold rj c = — Gi of the single mode 
solution. Part (b) of the figure holds an enlargement of the 
small rj range in order to distinguish properly the different 
curves in the lower half plane of part (a). Further parameters 
are p = 1.8 and 5=1. 

between the upper and lower solid lines is the portion 
of the 77-Gi-plane wherein one mode solutions are linear 
unstable. As predicted by Eq. I|29|) for p = 1.8 three 
mode solutions do not exist for modulation amplitudes 
greater than G\ — 0.3125. Before this critical value is 
reached the dotted and the solid lines are overlapping 
which means that the upper bifurcation of three mode 
solutions, i.e. the one occurring at large 77-values, has 
become supercritical as it has already been mentioned 
above. 

According to Eq. I jlOjl different amplitudes A and B 
generate dissimilar spatial patterns as illustrated in Fig.El 
for four distinct parameter sets. Fig. E^a) shows regular 
hexagons whereas in Fig.EJTb)-(d) distorted hexagons can 
be seen for various amplitude ratios A/B corresponding 
to the parameters used in Fig. 0(b) and (c). 
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FIG. 9: The spatial variation of one term contributing to 
Eq. 11011 is shown for hexagons in part (a) and for distorted 
hexagons in (b), (c) and (d). In part (b) a distorted hexagon 
is depicted for the parameters used in Fig. Ifjb) , namely Gi = 
0.05 and rj = -0.0492 with the amplitudes being A = 0.0551 
and B = 0.0124. In the remaining two parts the modulation 
amplitude is Gi = 0.11, whereas the control parameter is r\ = 
-0.099 in (c) (with A = 0.15 and B = 0.06) and 77 = -0.09 
in part (d) (with A = 0.198 and B = 0.116). 

C. Solutions of the forced amplitude equations for 

k m = k c and Gi =0, G 2 / 0. 

As in the previous case two amplitudes of the three 
mode solution coincide, i.e. Ai^ = Be 3,3 , whereas the 
third one is given by A\ = Ae tSl . The two remaining 
moduli A and B are determined by the following nonlin- 
ear equations 

(77- I A | 2 -2p I B | 2 ) A + G 2 + e m B* 2 = , (30a) 
(77 - (1 + p) I B | 2 -p I A | 2 ) B + e i@ A*B* = . (30b) 

If B = and Gi ^ the single mode solution is de- 
termined by a third order polynomial in A. In contrast 
to the case discussed in Sec. IV Bl thc bifurcation is now 
an imperfect one as can be seen in Fig. ^| Despite this 
imperfect bifurcation one may observe the following sim- 
ilarities to the previous paragraph wherein the 2 : 1 res- 
onance was studied: First of all the degeneracy of three 
mode solutions encountered for Gi — vanishes once 
again for finite values of the modulation amplitude Gi. 
Secondly, as for the 2 : 1 resonance, A corresponding 
to the largest amplitude of the three mode solution bi- 
furcates for small values of Gi from the top of the single 
mode solution, whereas the small amplitude B bifurcates 
from zero. 




-1.00 -0.50 0.00 0.50 fj -0.50 -0.10 0.30 0.70/1 



FIG. 10: The amplitudes of a three mode solution with A 
always corresponding to the upper branch and B to the lower 
one, as well as the amplitude of the imperfect bifurcating 
single mode solution (B = 0) are plotted as a function of the 
control parameter r\. Three different values of the modulation 
amplitude are taken into account: G2 = 0.05 (a), where (b) 
is an enlargement of (a) over a smaller 77-range, G2 = 0.1 (c) 
and G2 = 0.2 (d). Solid lines mark linear stable solutions and 
dashed lines unstable ones. Furthermore p = 1.7 and 5 = 1. 
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FIG. 11: In between the dotted line and the lower solid line 
(coinciding in the lower half-plane with the dashed line) three 
mode solutions exist. Along the dashed line, which overlaps 
on its lower branch with lower branch of the solid line, the 
values of the functionals of single and three mode solutions 
match. Thus within the area bounded by the dashed line 
three mode solutions are favored. Single mode solutions are 
linear stable in the complementary part of the region enclosed 
by the solid line. The nonlinear coefficient is assumed to be 
p = 1.7 and 5 = 1. 

With increasing values of Gi the bifurcation of the 
three mode solution changes once more at its lower bi- 
furcation point from a subcritical one in Fig. HOf a) to a 
supercritical one in Fig. HOf c) or in Fig. HOf d). The bi- 
furcation point located at large values of r\ is not shown 
in Fig. HOf a) and Fig. HOf c) but is subcritical as the one 
depicted in Fig. for a 2 : 1 resonance. Thirdly the 
hysteresis described by the three mode solution becomes 
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once again less pronounced with increasing values of Gi 
and turns out to be supercritical for large modulation 
amplitudes. Last but not least the three mode solution 
ceases to exist beyond some critical value of the modu- 
lation amplitude Gi as it was noticed as well for a 2 : 1 
resonance. 

The properties explained in the previous paragraph are 
recapitulated in Fig.^l where three mode solutions exist 
within the region bounded by the dotted lines, overlap- 
ping in the lower half-plane with the solid line. Along 
the dashed line the functionals T of the single and three 
mode solution match, meaning that the three mode so- 
lution is preferred in between the dashed lines coincid- 
ing once more with the solid line for small values of the 
control parameter r\. This overlapping of lines is due to 
the fact that the lower bifurcation point is supercriti- 
cal for small Gi and that the extension of the hysteresis 
shrinks with increasing modulation amplitudes. Through 
a linear stability analysis, similar to the one described in 
Sec. IIVCI one may conclude that single mode solutions 
are linear stable only outside of the area delimited by the 
solid line. 



VI. NUMERICAL RESULTS 

A major problem always encountered when dealing 
with amplitude equations is that their validity range 
around the threshold remains a priori unpredictable. 
Therefore we have determined the amplitude of the stripe 
solution by numerical simulations of the basic Eqs. (JTJ 
without considering modulations as a function of the rel- 
ative distance rj from threshold and compared these nu- 
merical results with the analytical solution A = y/rjpy in 
Fig- El Up to r\ ~ 0.2 we find a fairly good agreement be- 
tween analytical and numerical results, which gives a rea- 
sonable estimate of the quantitative validity range of the 
preceding results obtained in terms of amplitude equa- 
tions. 



VII. SUMMARY AND CONCLUSIONS 

The generic amplitude equations (|13[) are suggested for 
a forced two-dimensional system with a broken up-down 
symmetry and discussed throughout this work. Addi- 
tionally they are derived for the first time explicitly from 
a model for chemical reactions, in the present case the 
Lengyel-Epstein model. They are apt to describe how 
spatial periodic modulations interfere with the competi- 
tion between stripes and hexagonal patterns. Our inter- 
est was particularly focused on a 1 : 1 and a 2 : 1 reso- 
nance between the modulation wave number k m and the 
wave number k c of the respective pattern. Compared to 
the unmodulated case the bifurcation-scenarios of stripes 
and hexagonal patterns change strongly in the presence 
of modulations. It has been shown that three mode so- 
lutions, i.e. hexagons in the unmodulated case and dis- 
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FIG. 12: The amplitude of a stripe solution obtained by nu- 
merical simulations (circles) of the microscopic model given by 
Eqs. is compared with the analytical solution A = yrfpy. 
The parameters used are a = 16, a = 301, d = 1.07, c = 0.6 
which yield the nonlinear coefficient 7 = 0.288. For the sake 
of simplicity no modulations were taken into account. 



torted hexagons otherwise, are suppressed for large mod- 
ulation amplitudes. Furthermore evidence was given that 
depending on the strength of the modulation amplitude 
those three mode solutions bifurcate either sub- or su- 
percritically. The bifucation diagrams presented in this 
work may be confirmed by experiments on chemical re- 
actions with a spatially modulated illumination, where 
in contrast to Refs. |20l |2L| the ratio between the wave 
number of the illumination and the wave number of the 
Turing pattern should be chosen to be 2 : 1 or 1 : 1 . 

Instead of considering a stationary forcing the effect 
of traveling stripe forcing on the competition between 
hexagons and stripes has recently been investigated for a 
1 : 1 resonance in a one-dimensional system |34j . Given 
the coupled amplitude equations lj*P3|) the results estab- 
lished in Ref. |.3 1| may be directly interpreted within the 
scope of the microscopic Lengyel-Epstein model. The 
same allegation holds for the 1 : 1 resonance analyzed 
in a two-dimensional system in Ref. 35)]. Although the 
bifurcation-scenarios of the stationary 1 : 1 and 2 : 1 res- 
onance we studied in the present work share some simi- 
larities, there remains the open question how the bifur- 
cation diagrams presented in Refs. HH are modified 
if one has to deal with a 2 : 1 resonance instead of a 1 : 1 
resonance. 

We mainly focused on the spatially periodic forcing 
of Turing patterns. The Lengyel-Epstein model ex- 
hibits furthermore a Hopf-bifurcation 56] similar to those 
found in other chemical reaction schemes. Instead of a 
temporal resonant forcing of an oscillatory but spatially 
homogeneous chemical reaction |3ll l82j one may induce 
such an oscillatory reaction by a spatially periodic or a 
traveling stripe modulation. This leads to a revealing am- 
plitude equation, which has been recently deduced from 
the Lengyel-Epstein model and exhibits rich transition 
scenarios between solutions that are either spatially har- 
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monic or subharmonic with respect to the external mod- 
ulation Hal. 



This work results mainly from a summer project for 
students. 
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